Search for: All records

Abstract The inductive component of the magnetospheric electric field, which is associated with the temporal change of magnetic field, provides an additional means of local plasma energization and transport in addition to the electrostatic counterpart. This study examines the detailed response of the inner magnetosphere to inductive electric fields and the associated electric‐driven convection corresponding to different solar wind conditions. A novel modeling capability is employed to self‐consistently simulate the electromagnetic and plasma environment of the entire magnetospheric cavity. The explicit separation of the electric field by source (inductive vs. electrostatic) and subsequent implementation of inductive effects in the ring current model allow us to investigate, for the first time, the effect of the inductive electric field on the kinetics and evolution of the ring current system. The simulation results presented in this study demonstrate that the inductive component of the electric field is capable of providing an additional source for long‐lasting plasma drifts, which in turn significantly alter the trajectories of both thermal and energetic particles. Such changes in the plasma drift, which arise due to the inductive electric fields, further reshape the storm‐time ring current morphology and alter the degree of the ring current asymmetry, as well as the timing and the peak of the ion pressure. The total ion energy is increasing at a faster rate than the supply of energetic ions to the ring current, suggesting that the inductive electric field provides effective and accumulative local energization for the trapped ring current population without confining additional particles. 

Abstract Block-Adaptive-Tree Solar-wind Roe-type Upwind Scheme (BATSRUS), our state-of-the-art extended magnetohydrodynamic code, is the most used and one of the most resource-consuming models in the Space Weather Modeling Framework. It has always been our objective to improve its efficiency and speed with emerging techniques, such as GPU acceleration. To utilize the GPU nodes on modern supercomputers, we port BATSRUS to GPUs with the OpenACC API. Porting the code to a single GPU requires rewriting and optimizing the most used functionalities of the original code into a new solver, which accounts for around 1% of the entire program in length. To port it to multiple GPUs, we implement a new message-passing algorithm to support its unique block-adaptive grid feature. We conduct weak scaling tests on as many as 256 GPUs and find good performance. The program has 50%–60% parallel efficiency on up to 256 GPUs and up to 95% efficiency within a single node (four GPUs). Running large problems on more than one node has reduced efficiency due to hardware bottlenecks. We also demonstrate our ability to run representative magnetospheric simulations on GPUs. The performance for a single A100 GPU is about the same as 270 AMD “Rome” CPU cores (2.1 128-core nodes), and it runs 3.6 times faster than real time. The simulation can run 6.9 times faster than real time on four A100 GPUs. 

Abstract In a previous study, Huang et al. used the Alfvén Wave Solar atmosphere Model, one of the widely used solar wind models in the community, driven by ADAPT-GONG magnetograms to simulate the solar wind in the last solar cycle and found that the optimal Poynting flux parameter can be estimated from either the open field area or the average unsigned radial component of the magnetic field in the open field regions. It was also found that the average energy deposition rate (Poynting flux) in the open field regions is approximately constant. In the current study, we expand the previous work by using GONG magnetograms to simulate the solar wind for the same Carrington rotations and determine if the results are similar to the ones obtained with ADAPT-GONG magnetograms. Our results indicate that similar correlations can be obtained from the GONG maps. Moreover, we report that ADAPT-GONG magnetograms can consistently provide better comparisons with 1 au solar wind observations than GONG magnetograms, based on the best simulations selected by the minimum of the average curve distance for the solar wind speed and density. 

Abstract A potential field solution is widely used to extrapolate the coronal magnetic field above the Sun’s surface to a certain height. This model applies the current-free approximation and assumes that the magnetic field is entirely radial beyond the source surface height, which is defined as the radial distance from the center of the Sun. Even though the source surface is commonly specified at 2.5R_s(solar radii), previous studies have suggested that this value is not optimal in all cases. In this study, we propose a novel approach to specify the source surface height by comparing the areas of the open magnetic field regions from the potential field solution with predictions made by a magnetohydrodynamic model, in our case the Alfvén Wave Solar atmosphere Model. We find that the adjusted source surface height is significantly less than 2.5R_snear solar minimum and slightly larger than 2.5R_snear solar maximum. We also report that the adjusted source surface height can provide a better open flux agreement with the observations near the solar minimum, while the comparison near the solar maximum is slightly worse. 

Abstract The coronal heating problem has been a major challenge in solar physics, and a tremendous amount of effort has been made over the past several decades to solve it. In this paper, we aim at answering how the physical processes behind the Alfvén wave turbulent heating adopted in the Alfvén Wave Solar atmosphere Model (AWSoM) unfold in individual plasma loops in an active region (AR). We perform comprehensive investigations in a statistical manner on the wave dissipation and reflection, temperature distribution, heating scaling laws, and energy balance along the loops, providing in-depth insights into the energy allocation in the lower solar atmosphere. We demonstrate that our 3D global model with a physics-based phenomenological formulation for the Alfvén wave turbulent heating yields a heating rate exponentially decreasing from loop footpoints to top, which had been empirically assumed in the past literature. A detailed differential emission measure (DEM) analysis of the AR is also performed, and the simulation compares favorably with DEM curves obtained from Hinode/Extreme-ultraviolet Imaging Spectrometer observations. This is the first work to examine the detailed AR energetics of our AWSoM model with high numerical resolution and further demonstrates the capabilities of low-frequency Alfvén wave turbulent heating in producing realistic plasma properties and energetics in an AR. 

Abstract The vast size of the Sun’s heliosphere, combined with sparse spacecraft measurements over that large domain, makes numerical modeling a critical tool to predict solar wind conditions where there are no measurements. This study models the solar wind propagation in 2D using the BATSRUS MHD solver to form the MSWIM2D data set of solar wind in the outer heliosphere. Representing the solar wind from 1 to 75 au in the ecliptic plane, a continuous model run from 1995–present has been performed. The results are available for free athttp://csem.engin.umich.edu/mswim2d/. The web interface extracts output at desired locations and times. In addition to solar wind ions, the model includes neutrals coming from the interstellar medium to reproduce the slowing of the solar wind in the outer heliosphere and to extend the utility of the model to larger radial distances. The inclusion of neutral hydrogen is critical to recreating the solar wind accurately outside of ∼4 au. The inner boundary is filled by interpolating and time-shifting in situ observations from L1 and STEREO spacecraft when available. Using multiple spacecraft provides a more accurate boundary condition than a single spacecraft with time shifting alone. Validations of MSWIM2D are performed using MAVEN and New Horizons observations. The results demonstrate the efficacy of this model to propagate the solar wind to large distances and obtain practical, useful solar wind predictions. For example, the rms error of solar wind speed prediction at Mars is only 66 km s⁻¹and at Pluto is a mere 25 km s⁻¹. 

Abstract We describe our first attempt to systematically simulate the solar wind during different phases of the last solar cycle with the Alfvén Wave Solar atmosphere Model (AWSoM) developed at the University of Michigan. Key to this study is the determination of the optimal values of one of the most important input parameters of the model, the Poynting flux parameter, which prescribes the energy flux passing through the chromospheric boundary of the model in the form of Alfvén wave turbulence. It is found that the optimal value of the Poynting flux parameter is correlated with the area of the open magnetic field regions with the Spearman’s correlation coefficient of 0.96 and anticorrelated with the average unsigned radial component of the magnetic field with the Spearman’s correlation coefficient of −0.91. Moreover, the Poynting flux in the open field regions is approximately constant in the last solar cycle, which needs to be validated with observations and can shed light on how Alfvén wave turbulence accelerates the solar wind during different phases of the solar cycle. Our results can also be used to set the Poynting flux parameter for real-time solar wind simulations with AWSoM. 

Abstract Using a two‐way coupled magnetohydrodynamics with embedded kinetic physics model, we perform a substorm event simulation to study electron velocity distribution functions (VDFs) evolution associated with Bursty Bulk Flows (BBFs). The substorm was observed by Magnetospheric Multiscale satellite on 16 May 2017. The simulated BBF macroscopic characteristics and electron VDFs agree well with observations. The VDFs from the BBF tail to its dipolarization front (DF) during its earthward propagation are revealed and they show clear energization and heating. The electron pitch angle distributions (PADs) at the DF are also tracked, which show interesting energy dependent features. Lower energy electrons develop a “two‐hump” PAD while the higher energy ones show persist “pancake” distribution. Our study reveals for the first time the evolution of electron VDFs as a BBF moves earthward using a two‐way coupled global and kinetic model, and provides valuable contextual understanding for the interpretation of satellite observations.